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Abstract 


We investigate the cosmological attractor of the minimally coupled, self-interacting phan- 
tom field with a positive energy density but negative pressure. It is proved that the phantom 
cosmology is rigid in the sense that there exists a unique attractor solution. We plot the tra- 
jectories in the phase space numerically for the phantom field with three typical potentials. 
Phase portraits indicate that an initial kinetic term decays rapidly and the trajectories 
reach the unique attractor curve. We find that the curve corresponds to the slow-climb 
solution. 
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A scalar field with negative kinetic energy is proposed to explain the accelerated ex- 
pansion of present universe as dark energy [I]. This form of dark energy with the state 
equation parameter w < —1, dubbed phantom energy violates the dominant-energy condi- 
tion. It was shown that this model is consistent with both recent observations and classical 
tests of cosmology, in some cases providing a better fit than the more familiar models with 
w > —1. In spite of the fact that the field theory of phantom fields encounters the problem 
of stability which one could try to bypass by assuming them to be effective fields B] B], 
it is nevertheless interesting to study their cosmological implication. Recently, there are 
many relevant studies of phantom energy [4] and the primordial perturbation spectrum 
from various phantom inflation models [5]. 


The physical background for phantom type of matter with strongly negative pressure 
may be looked for in string theory [6]. Phantom field may also arise from a bulk viscous 
stress due to particle production [7] or in higher-order theories of gravity [8], Brans-Dicke 
and non-minimally coupled scalar field theories [9]. The cosmological models which allow for 
phantom matter appear naturally in the mirage cosmology of the braneworld scenario 
and in k-essence models [Ii]. 


In this letter we study the attractor behavior of phantom cosmology. Using the Hamilton- 
Jacobi formalism, we prove that there exists a unique attractor solution in the early universe 
containing a minimally coupled, self-interacting phantom field. We use an explicit numer- 
ical computation of the phase space trajectories for the phantom field with three typical 
potentials. Phase portraits indicate that the initial kinetic term decays rapidly and the 
trajectories reach the unique attractor curve. The attractor curve corresponds to the slow- 
climb solution. 


The action of the phantom field minimally coupled to gravity can be written as 
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where K? = 87Gy is the gravitational coupling, V(@) is the phantom potential and the 
metric signature (—, +, +, +) is employed. The climbing phantom in a spatially flat FRW 
cosmological model can be described by a fluid with a positive energy density p and a 
negative pressure P given by 
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p = —5¢°+V(9); (2) 


1 « 
= -58#-V(o), (3) 
which means that phantom energy violates the dominant energy condition, p + P < 0 and 
p > 0. The corresponding equation of state parameter is now given by 
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Since the phantom energy density is positive, Eq. (£) indicates that w < —1, and w — —1 
as the ratio ¢?/2V(¢) — 0. The evolution equation of the phantom field and the Friedmann 
constraint are 


+3H-— V'($) =0, (5 
m = £ -38+ vo). (6) 


Note that the sign of the potential force term is negative in Eq.(), which distinguishes 
the phantom field from the ordinary field and implies that the phantom field climbs up the 
potential. 


The Hamilton-Jacobi formulation is a powerful way of rewriting the equations of mo- 
tion, which allows an easier derivation of many inflation results. We concentrate here on 
the homogeneous situation as applied to spatially flat cosmologies, and demonstrate the 
stability of the phantom cosmology using the Hamilton-Jacobi formalism [2]. 


Differentiating Eq.(6) with respect to t and substituting in Eq.) gives 
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This allows us to write the Friedmann equation in the first-order form 
a2. 3 22 Kt 
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Eqs.(@) and (8) are the Hamilton-Jacobi equations. They allow us to consider H(¢), rather 
than V(¢), as the fundamental quantity to be specified. 


Suppose H (¢, p) is the solution of Eq. (8) which is uniquely determined once the initial 
conditions have been specified, where the parameter p is associated with each solution. The 
general solution to Eq.(@ can be expressed as 


donaa | [wmon (AP) | (9) 


where a; is the value at some initial point ¢;. We consider two solutions H (¢, p + Ap) and 
H(¢, p), where Ap < 1. By differentiating Eq.(8) with respect to p and using Eq. (Q), we 
obtain 


H(¢,p + Ap) — H(¢,p) x a™°(¢, p)Ap. (10) 


We find that any two solutions approach the attractor solution H(@) in an expanding 
universe. However, the attractor may not be the same for all values of the parameters. In 
order to check whether there exists a unique attractor, one defines the quantity [I3] 


r= E%r+Ap) l (11) 
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Eq.(10) shows that F — 0 as the universe expands, but the form of the attractor may vary 
if OF /0¢ changes sign. Note that F can go to zero for OF /0@ > 0 or for OF /0¢ < 0, but 
not for both. Thus, if OF'/0¢ = 0 for some value of the parameters, the attractor solution 
will not be unique for all values of the parameters. 
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which indicates that there exists a unique attractor because OF'/O¢ can not pass through 
Zero. 


To study an explicit numerical computation of the phase space trajectories, it is most 


convenient to rewrite the evolution Eqs.(5) and (6) for H and ¢ as a set of two first-order 
differential equations with two independent variables ¢ and @ 


w S, (13) 
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Let us consider three typical potentials for the phantom field: a quadratic potential 
V(o) = m?¢?/2, an exponential potential V(¢) = Voexp(—AK¢@) and a potential of the 
form V(¢) = Vo[l+cos(¢/f)]. We choose different initial conditions ġo and ġo in the range 
ba < 2V (ġo), and obtain the phase portraits (Fig.1, Fig.2 and Fig.3, respectively) in the 
(¢, $) phase plane. 


Power-law potential. In Figure 1, we see that there is a unique curve that attracts 
most of the trajectories in each branch of the two non-connected regions. An initial kinetic 
term decays rapidly. The phantom field climbs up the potential, which differs from the 
rolling-down behavior of the normal scalar field. The two attractor curves correspond to 
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cs 3k 

in the region ¢ > 0 and -mọ < ¢ < md, and 
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the slow-climb solutions: 


(15) 


in the region ¢ < 0 and mọ < ¢ < —md¢. In the model, the ratio of kinetic to potential 
energy of the phantom field tends to zero, so w — —1. The general features of the behavior 
of the trajectories do not change when we use a power-law potential V(¢) = A¢*. Although 
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Figure 1: Phase portrait for the phantom field with the potential V(¢) = m?¢?/2 in the 
(d, $) phase plane. 


Figure 2: Phase portrait for the phantom field with the potential V (¢) = Vo exp(—AK@) in 
the (¢,¢) phase plane. 


Figure 3: Phase portrait for the phantom field with the potential V(¢) = Vo[1 + cos(¢/f)] 
in the (¢,@) phase plane. 


for a > 2 the kinetic term is no more a constant value and increases as the universe expands, 
the ratio of kinetic to potential energy still tends to zero proportional to ¢~?. For a = 4 
the slow-climb regime becomes the regime of an exponential growth of the phantom field. 


Exponent potential. In Figure 2, the attractor curve corresponds to the slow-climb 


g=- [Penes (17) 


in the region —\/2V exp(—àÀkġ/2) < ġ < V2Voexp(—AKG/2). In the model, though the 
kinetic term grows exponentially the ratio of kinetic to potential energy of the phantom 
field is a constant value \?/6, so w = (A? + 6)/(A? — 6) < —1. Therefore, for a phantom 
field with an exponential potential V(¢) = Vo exp(—Ax«@), the scaling solution exists and is 
stable as long as A? < 6, while for \? > 6 the slow-climb approximation breaks down. This 


solution: 


regime is analogous to the scaling solution for a normal scalar field with an exponential 
potential [74]. 


Cosine potential. In Figure 3, the attractor curve corresponds to the slow-climb solution: 


b= a h + cos e) sin (18) 


in the region —/2Vo[1 + cos(o/f)]!/2 < ¢ < V2Vo[1 + cos(¢/f)]/?.. The phantom field 
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moves towards the top of the potential and then settles at o = 0. The friction term 3H db is 
so strong compared to the force term V’(¢) that the phantom field can not oscillate about 
the maximum of the potential, which differs from the regime for a gaussian potential [2]. In 
this case, the ratio of kinetic to potential energy of the phantom field tends to zero, so the 
state equation parameter w tends to —1. Thus the de-Sitter like solution is the late-time 
attractor of the model. 


In conclusion, we have investigated the attractor behavior of phantom cosmology. It 
is proved that there exists a unique attractor solution, hence the dynamical system of a 
single, minimally coupled phantom field is rigid irrespective of the nature of the potential. 
Moreover, we consider three different model: power-law potential, exponential potential and 
cosine potential, and plot the trajectories in the phase space numerically. Phase portraits 
indicate that an initial kinetic term decays rapidly and the trajectories reach the unique 
attractor curve. We find that the curve corresponds to the slow-climb solution. The work 
can be extended to the braneworld scenario. 
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